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Ginzburg-Landau theory is used to study the properties of single vortices and of the Abrikosov vortex 
lattice in a d x 2_ y 2 superconductor. For a single vortex, the s-wave order parameter has the expected 
four- lobe structure in a ring around the core and falls off like 1/r 2 at large distances. The topological 
structure of the s-wave order parameter consists of one counter-rotating unit vortex, centered at 
the core, surrounded by four symmetrically placed positive unit vortices. The Abrikosov lattice is 
shown to have a triangular structure close to T c and an oblique structure at lower temperatures. 
Comparison is made to recent neutron scattering data. 

PACS: 74.20.De, 74.60.Ec, 74.72.-h 
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Evidence continues to accumulate for the existence of 
nodes in the gap of high T c superconductors EL al- 
though the nature of the microscopic mechanism remains 
controversial. Thus there is good reason to study the 
phenomenology of superconductors with non-trivial gap 
structures, particularly their behavior in an applied mag- 
netic field that penetrates the superconductor, creating 
vortices which then play a dominant role in the transport 
properties of the system. We have previously considered 
0] a simple microscopic model of d-wave superconductiv- 
ity for electrons on a lattice and used the Bogoliubov-de 
Gennes equations to calculate the order parameter dis- 
tribution for a single vortex. The relevant Ginzburg- 
Landau (GL) free energy, which was first derived by 
Joynt ||, served to interpret our results. The GL the- 
ory involves both the d-wave order parameter and an 
induced s-wave order parameter which arises through a 
mixed gradient coupling. 
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Experience with conventional (s-wave) superconductors 
has demonstrated that virtually all of the phenomenolog- 
ical properties of superconductors can be derived from 
the appropriate GL theory. Here we present analogous 
results for what is arguably the simplest model of an 
unconventional superconductor with a non-trivial gap 
structure. First we integrate the GL equations for a sin- 
gle vortex numerically. The results, shown in Figure 1, 
exhibit a non-trivial internal topological structure for the 
s-wave component. Second we solve for the structure of 
the vortex lattice close to H C 2(T), by first minimizing 
the quadratic part of the free energy using a simple vari- 
ational wave function and then forming a periodic array 
of vortices from linear combinations of these functions. 
The results are compared to recent small angle neutron 
scattering (SANS) data j§. 

The free energy of a d x 2_ y 2 superconductor may be ex- 
pressed in terms of the two order parameters, s(r) and 
d(r), with appropriate symmetries, as follows || 

(3 3 \s\ 2 \d\ 2 +(3 4 {s* 2 d 2 +d* 2 s 2 ) 

U yS )*(Uyd) - (U X S)*(U X d) + C.C.] . (1) 



Here n = —ih\7 — e*A/c, and d is assumed to be the critical order parameter, i.e., we take a s = T — T s , ay = T — Td 
with T s <Td- It is also assumed that f3\, /3%, (3$, j s , 7^ and j v are all positive B. The parameters ji are related to the 
effective masses in the usual way. We use 7$ = h 2 /2m*, for i — s,d,v. The field equations for the order parameters 
are obtained by varying the free energy (hi) with respect to conjugate fields d* and s*, giving, 
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(n 2 -n 2 )s + 2/3 2 |d| 2 d + /3 3 |s| 2 d + 2/3 4 s 2 d* = 0, 
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-(II 2 - n 2 )d + 2(3 l \s 2 \s + /3 3 |d| 2 s + 2/? 4 cZ 2 s* = 0. 
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Equations (||) can be integrated numerically for boundary conditions which generate a single vortex at the origin. 
In doing this, we assume an extreme type-II limit, where the coupling to the vector potential can be ignored while 
considering the core structure of the isolated vortex line. 
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Ren et al. [El have previously shown that for a d-wave 
order parameter with the asymptotic form, 



d(r, 9) = d e ie 



(3) 



where do = y/—&d/2P2, the asymptotic form of the s- 
wave order parameter is: 



s(r,9) = 9i(r)e 



(4) 



where g\ir) and gi(r) fall off like 1/r 2 for large r. Fur- 
thermore, close to Td, g2 (r) ~ —3gx(r) and therefore the 
winding number far from the core is +3. This result com- 
bined with the result that close to the core the winding 
number is — 1 0] forces us to conclude that four addi- 
tional positive vortices must exist outside the core. We 
emphasize that this is a topological result and thus not 
sensitive to small modifications of the parameters. As is 
shown below, these vortices lie on the ±x and ±y axes. 
At lower temperatures a topological transition to a state 
with s-wave winding number —1 is in principle possible. 

Asymptotically the superconductor is not in a pure d- 
wave state, but rather in a state characterized by power 
law decay of the s-wave component. Only at the length 
scale given by the penetration depth is the pure d-wave 
state regained. 

We have studied the dependence of the maximum of 
the s-wave component on the GL parameters. Noting 
that both the d-wave and s-wave components rise over 
the same length scale given by £d, where = 7d/|ay|, 
allows us to give an order of magnitude estimate for the 
magnitude of the s-wave order parameter at the maxi- 
mum, 
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Our numerical results confirm that the constant of pro- 
portionality is of the order unity. Note that the temper- 
ature dependence of max(s) is (1 — T/Td) 3 / 2 . 

In Figure [l] we show the behavior of the s-wave am- 
plitude along the x-axis and along the diagonal, as ob- 
tained by numerical integration of Eqs.(j|). Moving out- 
ward from the center of the vortex both the s- and d-wave 
amplitudes increase over the same length scale £d- Mov- 
ing further out one enters a region where the relative 
phase tends to lock to value ±7r/2. The change in the 
relative phase takes place in narrow "domain walls" . Up 
to this distance the results are in perfect agreement with 
the ones obtained within Bogoliubov-de Gennes theory 
||. However, further out the situation changes. The 
domains of rapid variation of the relative phase vanish. 
Furthermore, the relative phase starts to wind in the op- 
posite direction. In Figure [l] this change manifests itself 
as a zero in the amplitude of the s-wave component. This 
zero is nothing but the core of one of 4 "extra" vortices. 
Identical vortices are found at all the four "wall ends"; 



this combined with the vortex with an opposite charge 
at the center gives the total required winding of +3. 

Next we turn to the problem of the structure of the 
vortex lattice in the vicinity of the upper critical field 
H C 2 where the amplitudes of the order parameters are 
small and it is sufficient to consider the linearized GL 
equations. It is easily seen that in the Landau gauge 
(A = yBx) these linearized field equations are satisfied 
by taking d(r) — e tky d(x), s(r) = e zky s(x). Then, ex- 
actly as in the one component case fij] , we are left with a 
one dimensional problem which can be stated as follows: 



{Ho + a d )d + Vs = Ed, 



Vd+(H + a s )s = Es, 
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where 7io = huj c (a^ a-\-\/2) and V = e v (?iuj c /2)(a f a^+aa) 
are expressed in terms of the usual raising and lowering 
operators, which can be written as a = [(x — Xk)/l + 
l(d/dx)]/y/2. Here I = y / hc/e*B is the magnetic length, 
Xfe = kl 2 and oj c — (e*B/mc). In writing Eqs. (Q), we 
have assumed, for simplicity, that m* d = m* = m , i.e. 
that ^d — Is, and we have set e v = ^vHs = m t/ m l- By 
including the right hand side of Eqs.(|6j) we are consider- 
ing a slightly more general problem: E = corresponds 
to the physical solution for B = H C 2 (T) , and solutions for 
E < will be useful later when we consider the stability 
of various vortex lattice structures. 

In contrast to the one component case, the linearized 
Eqs.(0) have no obvious exact solutions. This is due 
to the coupling term V whose origin traces back to the 
mixed gradient term in the free energy (0). In what fol- 
lows we construct a simple variational solution, which is 
likely to capture all the essential physics of the problem. 
To this end we define H.^ = Ho ± V, and ip^ = d ± s, in 
terms of which we can write the set of equations (^|) as 



H+ +T -T* -AT/2 




where we have defined T* = (T d + T s )/2, AT = T d -T s . 
A nice feature of this representation is that for AT = 
the equations for ip + and <p~ decouple, each becoming 
a simple harmonic oscillator problem. Motivated by this 
fact, we consider a variational solution to the full problem 
in terms of normalized ground state harmonic oscillator 
wave- functions , 



z (x-x k ) 2 /2 
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The variational parameters a + and cr_ will be determined 
by minimization of the eigenvalue (E) . If o~ + = a cos i? 
and <7_ = (7sini9, the resulting expression for (E) can be 
minimized with respect to a 2 analytically. This leads to 
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where x = tan$. The last equation must be minimized 
numerically with respect to x, and is governed by the 
parameter A = tku c /AT. In the low field limit, set by 
Acl, (7+ ~ <r_ ~ ^/muj c /h, while in the high field limit, 
A > 1, cr± w ^mJ c /ft[(l ± e„)/(l t e*)] 1/4 - It follows 
that at least intermediate values of A are required for 
appreciable effects from s-d mixing to occur. Otherwise 
the s component effectively vanishes. 

Solutions to Eq.(||) with (E) = give the dependence 
of the upper critical field H c2 on the temperature. When- 
ever a finite admixture of the s component is present, a 
characteristic upward curvature is found near the critical 
temperature in H C 2(T). Such curvature has been exper- 
imentally observed in both La-Sr-Cu-0 and Y-Ba-Cu-0 
compounds || , and has also been interpreted as a result 
of s-d mixing by Joynt 

We next construct a vortex lattice. Consider a periodic 
solution of the form 



Xd/s(r) = ^2c n e l 



™[<p+(x)±<p;(x)], 



(9) 



where f^(x) are defined by (Q) and k = qn, (n integer), 
which gives periodicity in y with period L y = 2n/q. So- 
lution (^) will also be periodic in x provided that the 
constants c„ satisfy the condition c n+ N = c n for some 
integer N. In what follows we consider only the case of 
N = 2 so that C2n — cq and C2 n +\ = c\- The period 
in the x direction is L x — 2l 2 q, and it also follows that 
BL x L y = 2(hc/e*) — 2$o, i.e. there are two flux quanta 
per unit cell. The resulting lattice may be thought of 
as centered rectangular with two quanta per unit cell 
or, equivalently, as an oblique lattice with lattice vec- 
tors of equal length and one flux quantum per unit cell. 
The parameter q controls the shape of the vortex lattice, 
and it is customary to describe this shape by the ratio, 
R = L x /L y = (I 2 /ir)q 2 . R = 1 corresponds to the square, 
while R = a/3 corresponds to the triangular vortex lat- 
tice. The restriction to centered rectangular lattices is 
made primarily for computational convenience. However 
it is compatible with recent experiments on YBa2Cu3 07 
which show evidence for an oblique vortex lattice with 
nearly equal lattice constants Q . 

At B ss H c2 (T) all solutions of the form (||) are degen- 
erate. It is the fourth order terms in the free energy that 
lift this degeneracy below H C 2 and determine the vortex 
lattice configuration. Consider the solution to the full 
nonlinear problem of the form (s,rf) = c(x s ,Xd) where 
c is an arbitrary constant. Then the total integrated 
free energy of the system becomes F — |c| 2 (/2[Xs, Xd]) + 



volume of the system, and /2 and fa stand for quadratic 
and quartic parts of the free energy density (|l|) respec- 
tively. One can now minimize the total free energy F 
with respect to |c| 2 to obtain 
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m = 0A ' 
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l (fa[Xs,Xd}}, where 



means integration over the 



where (3a is the generalization of the usual Abrikosov 
parameter (t]]. 

We have studied the dependence of (3a on R in vari- 
ous regions of the parameter space. The qualitative fea- 
tures of this dependence are most strongly influenced by 
the parameter e v . In particular, for values e v close to 
0, (3a(R) is minimized by R m in = \/3, i.e. the triangu- 
lar lattice is stabilized. This is to be expected, since in 
this limit o+ ~ <r_, so the s-component is suppressed 
and the usual one component solution is found. 

The consistency of our solution is confirmed in this limit 
where we recover the correct value of the Abrikosov ratio 
@a(V3) = 1.1596 as quoted by Kleiner et al. How- 
ever as e v is increased, the minimum of (3a moves toward 
smaller values of R m in < \/3 signaling that an oblique 
vortex lattice is preferred. Finally, at some value of e v 
(which depends on other parameters) the minimum of 
Pa reaches Rmin = 1, characteristic of the square vortex 
lattice. Further increase of e v has no effect on the lattice 
structure which remains square. The typical dependence 
of (3a on R is displayed in Figure g. 

For all different parameter combinations studied, the 
minimum R m in varies continuously with e v . Varying 
other parameters changes somewhat the quantitative be- 
havior of this dependence, but the essential qualitative 
features described above remain intact. For example, in- 
creasing the parameter (3± makes the transition from tri- 
angular to square lattice somewhat sharper (see inset of 
Figure ||). One would also expect the inclusion of three 
dimensional screening effects (i.e. replacing B = H with 
B = H + (h s ), where (h s ) is calculated self-consistently) 
to change the results quantitatively, but not qualitatively. 

In practice this continuous dependence would mean 
that in a ei-wave superconductor one would expect to ob- 
serve a general oblique vortex lattice, unless the material 
is in one of the limiting regimes where e v is very small or 
very large. Such an oblique lattice structure has in fact 
been recently observed by SANS on YBa 2 Cu30y single 
crystals in strong magnetic fields parallel to the c axis 
by Keimer et al. Q. These authors reported an oblique 
lattice structure with nearly equal lattice constants and 
an angle of <f> = 73° between primitive vectors. Our phe- 
nomenological theory is consistent with any angle <f> in 
the interval (60°, 90°), including that found experimen- 
tally. We display an example of a general oblique vortex 
lattice obtained by explicitly evaluating amplitudes of s 
and d components of the order parameter from Eqs. (|^) 
in Figure ||. Comparison of Figures 0a and ||b reveals that 
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the non-trivial topological structure of the s component 
persists even in this high field regime. 

Keimer et al. further report that one principal axis of 
the oblique unit cell is always found to coincide with the 
(110) or (110) direction of the YBa 2 Cu 3 07 crystal. This 
is at variance with our results, since we find that one lat- 
tice vector of the larger rectangular cell is oriented along 
(100) or (010), even in the presence of a small orthorhom- 
bic distortion. However, we note that the energy cost of 
rotating the vortex lattice is small compared to the en- 
ergy needed to deform the lattice. It is thus possible 
that (110) twin boundaries, where the order parameter 
is weakened, bind lines of vortices and hence orient one 
of the oblique lattice vectors along (110) as is found ex- 
perimentally. 
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FIG. 3. Contour plot of the amplitudes of (a) d component 
and (b) s component of the order parameter. GL parameters 
are the same as in Figure ^ with e v — 0.45, resulting in an 
oblique vortex lattice with Rmin = 1-29 and the angle between 
primitive vectors <j> = 76° . The lightest regions correspond to 
the largest amplitudes. 
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FIG. 1. Amplitude of the s-wave component along the 
re-axis (solid line) and along the diagonal x — y (dotted line) 
normalized to the bulk value do- The parameters used are: 
7s = 7d = lv, ci s = 10|a d |, 0i = fa = 0, and 04, = 0.5/3 2 - The 
inset shows schematically the positions of the s-wave vortices 
and their relative windings. 



FIG. 2. Abrikosov ratio /3a as a function of the lattice 
geometry factor R — L x /L y for different values of e v and 
T s = 0.5T d , T = 0.75T d , ft = ft = ft = ft = 1, and 
B = 0.8H C 2- The inset shows the dependence of the mini- 
mum Rmin on the parameter e„ for different values of ft. 
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